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It has been hypothesized recently that core collapse supernovae are triggered by mildly 
relativistic jets following observations of radio properties of these explosions. Association 
of a jet, similar to a gamma-ray burst jet but only slower, allows shock acceleration of 
particles to high energy and non-thermal neutrino emission from a supernova. Detection 
of these high energy neutrinos in upcoming kilometer scale Cherenkov detectors may be 
the only direct way to probe inside these astrophysical phenomena as electromagnetic 
radiation is thermal and contains little information. Calculation of high energy neutrino 
signal from a simple and slow jet model buried inside the pre-supernova star is reviewed 
here. The detection prospect of these neutrinos in water or ice detector is also discussed 
in this brief review. Jetted core collapse supernovae in nearby galaxies may provide the 
strongest high energy neutrino signal from point sources. 
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1. Introduction 

Core collapse of massive stars which lead to supernovae (SNe) of type Ib,c and II are 
in some cases associated with long duration (> 2-10 3 s) gamma-ray bursts (GRBs), 
as evidenced by observed correlations of GRB 980425/SN 1998bw, GRB 021211/SN 
20021t, GRB 030329/SN 2003dh and GRB 0131203/SN 20031wP A relativistic jet 
with bulk Lorentz factor Tf, > 100, powered by a black hole and an accretion disc 
which form after the core collapse in the most likely scenario, is believed to lead 
to the GRB event Observational evidence of only a small fraction of detected 
SNe associated with GRBs hints that the frequency of highly relativistic jets in 
core collapse SNe is at best 1 in 1000, roughly the ratio of GRB to SN rates. ^ 
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However, a significantly larger fraction (< 10% of type Ib/c ratePEj) of S Ne (also 

called hypernovae) may have mildly relativistic jets associated with 

One or more of the following observations support the jetted SN hypothesis: 

• High expansion velocity (30-40 xlOOO km/s) first o bse rved in SN 1998bwPI 

• Radio afterglow not associated with 7-ray emission.^^ 

• Asy mmetr ic explosion supported by polarimetry observations of SN type 

ib/c.™i 

Numerical simulations of core collapse SNe, carried out over the last three 
decades have failed to produce a success ful explosion by a prompt shock wave cre- 
ated due to the collapse of its iron core^^ The deposition of bulk kinetic energy 
in a jet for m into the stellar envelope may help disrupt and blow it up making the 
SN possiblePEI T he presence of a jet is also conducive to shock acceleration of 
particles. In case of a GRB, internal shocks of plasma material along the jet accel- 
erate protons and electrons which radiate observed 7-raysP^ High energy protons 
may escape a s co smic rays and/or produce 100 TeV neutrinos by interacting with 
7-rays in situ\^ While the GRB jet is making its way out of the collapsing stellar 
progenitor it is expected to produce 10 TeV neutrio precursor burst (lSlISI These 
neutrinos are emitted even in the cases when the jets do not manage to burrow 
through the stellar envelope and choke inside without producing observable 7-rays. 
The jets in core collapse SNe or hypernovae which is the topic of this review are 
slow with Tft ~ few and choke inside the stellar envelope !^ Ne utrinos produce from 
such jets are typically of a hundred GeV to TeV energy! u |1 

As opposed to 10 MeV thermal neutrinos produced by the core collapse SN 
shocks which have been detected from SN 1987A in our own galaxyE^I a high en- 
ergy neutrinos from the jets may be detected from a longer distance because of an 
increasing detection prospect with n eutr ino energy. Kil om eter scale ice and water 
Cherenkov detectors such as IceCube^ and ANTARES^l which are currently be- 
ing built in Antarctica and in the Mediterranean will have an excellent chance to 
detect these neutrinos from SNe within the nearest 20 Mpc. 

The organization of this brief review is as follows: In Sec. [21 a basic core collapse 
SN picture is outlined and a particular slow jet model in Sec. El Shock acceleration 
and the maximum energy reachable by protons are discussed in Sec. Neutrino 
flux on Earth from a point source and diffuse sources is calculated in Sec. [5] and 
their detection prospects in Sec. Conclusions are given in Sec. 

2. Core Collapse Scenario 

Nuclear fusion reactions, similar to the ones which take place in our sun, constantly 
enrich the interior of a star forming an iron core as the end product. Burning up all 
fusion materials causes hydrodynamic instability due to lack of radiation pressure 

a see Ref. 1231 for a review of current experimental status to detect these low energy neutrinos. 
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from inside the star. The immense gravitational pressure of the stellar envelope 
and/or overlying material causes the core of stars with mass > 8Mq to collapse 
at this point. The density of the compressed core material reaches a few times the 
nuclear density and a rising temperature helps iron dissociate into nucleons and 
alpha particles. 

Infall of stellar material onto the core produces ~ 10 MeV neutrinos by the 
process of electron capture on protons [e~p — > v e n). The density of neutrons in 
the core exceeds that of protons in this process, called neutronization. Initially the 
neutrinos are trapped within a radius called neutrino sphere because of a density 

> 10 12 g-cm~ 3 . For progenitors of mass < 28M , the increasing degeneracy pressure 
of the neutrons leads to a rebound, which sends a Shockwave through the core. 
While traversing through the core, the Shockwave heats up material, dissociates 
more iron atoms and releases trapped v e from the neutrinosphere. Neutrinos carry 
away < 10 52 erg of energy or roughly 1% of the total gravitational binding energy 
in this bursting phase which lasts for a few milliseconds. The Shockwave, however, 
does not reach the envelope to drive it away because of heavy energy loss and the 
star fails to explode into a supernova. 

The mechanisms envisaged to produce a successful supernova explosion, such as 
observed in nature, may be divided into two main categories despite many uncer- 
tainties such as the mass loss rate of the pre-supernova star and neutrino transport 
in the core, to name a few. The first is a revived shock model, for stars initially less 
massive than ~ 28M Q , where the core collapses to make a neutron star. In this 
case, the above-mentioned stalled supernova shock is re-energized by neutrino ab- 
sorption on nucleons outside the stellar core (y e p — > e + n; v e n — > e~p), re-energizing 
them. The shock wave then reaches the envelope and expels it away. After the su- 
pernova explosion, the stellar core cools down in next 10's of seconds by emitting 
~ 3 x 10 53 erg of energy in neutrinos of all flavors created by lepton pair annihila- 
tion (e + e~ — » v^.tV^.t), neutrino pair annihilation (v e v e — > v^^v^^) and nucleon 
bremsstrahlung (NN — > NNv^^v^^). These neutrinos are thermal with a typical 
energy ~ 10 MeV. A neutron star is left over after the stellar core cools down, 
following ejection of the envelope and outer core. 

In the second scenario, which is the relevant case for the current review, a star 
initially more massive than ~ 28M Q undergoes core collapse, (a) If the initial mass 
is 28M Q < M* < 40M Q , the core collapses initially to a neutron star, but after fall- 
back of additional core gas which did not reach ejection velocity, it collapses further 
to a black hole (BH) of mass in excess of the Chandrasekhar mass (the neutron 
degeneracy pressure not being sufficient to counteract gravity for this mass), (b) 
If the initial mass M* > 40 Mq, the core collapses directly to a black hole of mass 

> 3 — 5Mq (see Ref. l26l e.g.). Thermal neutrinos of ~ 10 — 30 MeV are also produced 
in the neutronized core, before it falls into the black hole. This neutrino luminosity, 
of order several solar luminosities, may or may not be able to eject (via absorption 
in re-energization) the outer envelope, which is needed in order for it to appear 
as supernova detectable by its photon emission. As in the case of stars of mass 



February 5, 2008 7:55 WSPC/INSTRUCTION FILE snc 



4 Razzaque, Meszdros, & Waxman 

< 28M &1 numerical simulations have not yet been able to prove whether stars in 
this mass range eject their envelope; and, observationally, it is unclear whether any 
observed supernova can be ascribed to progenitors in this mass range. However, 
stars more massive than ~ 28 Mq are certainly observed, and from well understood 
physics, they must core-collapse to a black hole. This is the gist of what happens 
to stars of mass > 28M Q if the stellar core is rotating slowly. A black hole, and 
possibly a temporary small accretion disk is formed after core collapse, which may 
not greatly affect the symmetry of the collapse and/or envelope ejection. 

The situation is thought to be drastically different for stars in the range M* > 
28M© whose core is fast-rotating P3 In particular, for a core angular momentum 
in the range (3-20) xlO 16 cm 2 s _1 for this model, the black hole and disk accretion 
can serve as an energy source for powering a long GRB, and may be able to eject 
a stellar envelope with kinetic energies possibly 10-100 times higher (at least in 
an isotropic-equivalent sense) than the energy of typical type Ib/c or type II SNe. 
Fast core rotation also helps in forming low density channels along the rotation 
axis, by centrifugal evacuation. Part of the material accreting from the disc onto 
the black hole can then be ejected as narrow jets along the axes, collimated by 
the gas pressure of the envelope. The jets are powered by neutrino annihilation or 
magnetohydrodynamic stresses. The pressure and energy deposition from the jet 
helps eject the star's envelope and thus the SN explosion happens. In the case of 
a highly relativistic jet which manages to break through the stellar envelope, the 
GRB event takes place outside the star in an optically thin environment. A slow 
jet, which is modelled in the next section, may never break through, but the ejected 
material can give rise to an "orphan" radio afterglow which is not associated with 
7-ray emission. However, non-thermal high energy (~ TeV) neutrinos are produced 
ubiquitously in both cases. 

3. A Slow Jet Model 

The pre-supernova star, after losing its outermost envelope, is typically a Wolf- 
Rayet star of radius ~ 10 11 cm in the case of a type lb or larger in the case of 
a type II supernova. The mildly relativistic SN jet may be modelled inside the 
pre-supernova star (see Fig. with a bulk Lorentz factor Tb ~ 10°" 5 r& 0.5 an d a 
total jet kinetic energy Ej ~ 10 51 ' 5 ^5i.5 ergs which is ~ 1% of the total energy 
released in the SN explosion. Because of its relativistic motion, the jet is most likely 
beamed with an opening half angle 6j ~ 1 /Tb which is much wider compared to a 
GRB jet of Tb > 100. Assuming a duration tj ~ lOt^i s, the isotropic equivalent 
kinetic luminosity of the jet is Lkin — Z^lEj/tj. With a jet variability time scale 
t v ~ O.li^-i s at the base, internal shocks between plasma materials moving along 
the jet occur at a radius rj ~ 1T\ct v w 10 10 8 r 2 5 t„,-i cm, which is below the 
stellar surface. 

Internal shocks convert a fraction e e ~ 0.1e e ,-i of the bulk kinetic energy (Ej) 
into random electron motion analogous to the GRB fireball models. In the case of 
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Fig. 1. Cartoon of a mildly relativistic jet buried inside the envelope of a collapsing star. The jet 
may deposit its kinetic energy and help blow up the envelope but may not get out itself. 



GRBs, these relativistic electrons would emit synchrotron photons in an optically 
thin environment, which are observed as 7-rays on Earth. The density of these 
electrons and baryons (since they are coupled to the baryons) in the jet may be 
estimated as 

■•' " ■■' - ^ - E] ~ 1Q 20.5 £51.5 cm _ 3j (|) 



4irr 2 r 2 m p c 3 2Trr 2 j m p cH j ' ' 11,0.6^.1*5 

b 



in the comoving jet frame. The opacity to Thomson scattering by photons 

t , VVzKZl « 1 ^_^ (2) 



in the comoving frame is then very high. 

The magnetic fields in the jet, built up by turbulent motions in the shock region, 
have an energy characterized as a fraction eb ~ O.Isb _i of the total jet energy, 
B 2 /8tt ~ £sikin/47rr|r^c. The corresponding magnetic field strength in the jet is 
given by 



B ^l^f~ w9 ( r*'lT V /2 «- (3) 



Electrons and protons are expected to be accelerated to high energies in the internal 
shocks, via the Fermi mechanism. The electrons cool down rapidly by synchrotron 
radiation in the presence of the magnetic field in Eq. ©. However, due to the large 
Thomson optical depth in Eq. J5J, these photons thermalize and the corresponding 
black-body temperature is 



b Hereafter the variables in the comoving and in the laboratory frame or local rest frame are 
denoted with and without a prime respectively. 
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The volume number density of these thermal photons may be roughly calculated as 

/ 3 3 \ !/ 4 

m > ~ £e Ej ^m24.8/ £ e, -1-^51. 5 \ „ TO -3 



„ ^ ~ 10 24 « . - - cm \ (5) 

It may be noted that photons from the shocked stellar plasma do not diffuse into 
the jet due to the high optical depth r^ h in Eq. @, and the number density in Eq. 
® is roughly constant. 



4. Proton Acceleration and Cooling Processes 

The shock acceleration time for a proton of energy E' p is proportional to its Larmor's 
radius and may be estimated as 

where k ~ 10«i and the magnetic field in Eq. (J^J was used. The maximum proton 
energy is limited by requiring this time not to exceed the dynamic time scale for the 
shock to cross plasma material: t' dyn ~ t v Tb, or any other possible proton cooling 
process time scale which we discuss next. 



4.1. Electromagnetic cooling channels 
4.1.1. Synchrotron and inverse Compton 

The cooling time scale for protons by synchrotron radiation in the same magnetic 
field which is responsible for its acceleration is given by 



^m p c 3 / E' p \ 1 / r,!..,,/,. : / 



t > El w 3 8 | ) I zM^lSzl s ( 7 ) 

~ a Th PmlE> p B'2 \GcV J { e B ^E 51 . 5 J S ' [ '> 

with (3 — v / c ~ 1. Inverse Compton (IC) scattering of thermal electron synchrotron 
photons is another cooling channel for high energy protons. The IC cooling time 
scale in the Thomson and Klein-Nishina (KN) regimes, valid for E' p much less or 
greater than m^/E^ w 10 5 - 3 (r^ 5 tj,itl _ 1 /E^i^e e -i) 1 / i GeV respectively, as 



3 m y oa f K Y ( T UJi,A-i \ 

'"■ n ' " -^—^ ^{Gev) { e e ^E 51 . 5 ) 

1/2 



V 7 7 



51.5 



Here we used the thermal photons with peak energy and density given in Eqs. (QJ 
and respectively as targets. 
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4.1.2. Bethe-Heitler 

Because of a high density of thermal photons in the SN jet in Eq. (j^J, protons 
may produce e + e~ pairs by interacting with them, a process known as Bcthe- 
Heitler (BH). The cross-section for BH interaction: jry — > pe + e~ is given by ctbh = 
ar 2 {{2S/9)\n[2E' p E' 1 /(rn p m e c iL )] - 106/9}. The logarithmic increase of the cross- 
section with incident proton energy implies that this is a very efficient cooling 
mechanism for high energy protons. The e + e~ pairs are produced at rest in the 
center of mass (cm.) frame of the collision and acquire an energy m e c 2 j' c m each in 
the comoving frame. Here i c m = (E' p +E! y )/(m 2 p c' i + 2E' p E! y )^ 2 is the Lorentz boost 
factor of the cm. in the comoving frame. The energy lost by the proton in each BH 
interaction is thus the energy of the created pairs in the comoving frame, and is 
given by AE' p = 2m e c 2 j' c m . The energy loss rate of the proton is proportional to 
the BH scattering rate as given by dE' p /dt' BH — n' cobhA-E' , and the corresponding 
proton cooling time is 

,, K E' p (m 2 p c* + 2E' p E>y/ 2 

BH dE p /dt' BB 2^a BH ™ e c 3 (^+-B;)' {) 

in the comoving frame. 



4.2. Hadronic cooling channels 

Photomeson (jry) and proton-proton (pp) interactions which are responsible for 
producing high energy neutrinos may also serve as a cooling mechanism for the 
shock accelerated protons. The p"f at the A + resonance and the average pp cross- 
sections are a vl = 5 x 10~ 28 cm 2 and a pp « 5 x 10 -26 cm 2 respectively. The 
corresponding optical depths, given by 

_ PyyT^rj 78 / eg_i^|i. 5 \ f 
V/ - 7T ~ tu -g -3—2 and 
1 b 6,0.5^,1 L v,-i ) 

I a PP n 'p r 3 ^ 1n 5.5 ( E 51.5 ^ nn s 

pp = "TT - r 3 t ,t 1 ' ( 0) 

are very high. The threshold proton energy for A + production against the target 
thermal photons of energy E' in Eq. (3J is 

E> A+ = ^! „ 10 - ] GeV. (11) 



T 



51.5 




Adopting the energy loss by a proton AE' p w 0.2_Ep and 0.8E' p respectively for 
each j»7 and pp interaction, the hadronic cooling time scales are 

rpl /r 12 +3 j-6 \ 1 / 4 

t > - io-T-3 W!d s 

" ca p ^AE' p \ 4,-iE$t. 5 
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ca pp n' p AE p 




10 



-5.6 




(12) 



using Eqs. JSJ and (JTJ. The photomeson cooling time scale t' pl above is roughly valid 
at E' = E' » + and one needs to use the thermal photon spectrum to calculate it 



4.3. Maximum proton energy 

The shock acceleration time in Eq. © for protons in the SN jet (solid line) and 
the different cooling time scales are plotted in Fig. [3] as functions of the comoving 
proton energy. The dashed lines are hadronic cooling time scales in Eq. I|12(l c and 
the dot-dashed lines are electromagnetic cooling time scales in Eqs. J7J, (JSJ and 
©. Note that the hadronic cooling time (t' pp ), the BH cooling time (t' BE ) and the 
synchrotron cooling time (t' syn ) are first longer and then shorter than the maxi- 
mum proton acceleration time (^ c ). The hadronic (i' ) and IC scattering (t' IC ) are 
not efficient cooling mechanisms for protons. The maximum proton energy can be 
roughly estimated, by equating the t' to t' accl from Eqs. © and Q as 



since t' BH w t' syn w t' acc at this energy. 

5. Neutrino Production and Flux on Earth 

Shock accelerated protons in the SN jet can produce non-thermal neutrinos by 
photomeson (jry) interactions with thermal synchrotron photons and/or by proton- 
proton (pp) interactions with cold protons present in the shock region. As shown in 
Fig. [21 the pj process is dominant in the energy range E' p w 10 4 - 10 52 GeV and 
the pp process is dominant at other energies. 

In the case of pj interactions at the A + resonance, neutrinos are produced 
from charged pion decay as pj — > A + — > mr + — > p + v^ — > e + v e i>^v^. The 

pp interactions also produce charged pions (tt^) and kaons [K ) and their decay 
modes are the same as above with 100% and 63% branching ratios respectively. d 
The total pion (kaon) multiplicity in eachpp interaction is ~ 1 (<~ 0.1) in the energy 
range considered hereP^The energy of the shock accelerated protons in the SN jet 
is expected to be distributed as oc 1/ E p , following the standard shock acceleration 
models. Charged mesons, produced by pp andp7 interactions, are expected to follow 
the proton spectrum with ~ 20% of the proton energy for each pion or kaon. 

C A thermal photon spectrum is used to calculate p7 cooling with a delta-function approximation. 
d 7r~ and K~ decay modes are charge conjugate of 7r+ and decay modes. 



txb n/ p — rj ^ + ciiiu une neei 

for different proton energies. 




(13) 
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Proton energy E' p (GeV) 

Fig. 2. Proton cooling time scales for different hadronic (dashed lines) and electromagnetic (dot- 
dashed lines) processes: photomeson (t' pl ), proton-proton (t' p ), synchrotron radiation (t' Byn ), IC 
scattering (tj c ) and BH (ig H ); in the comoving frame, as functions of energy. The shock accel- 
eration time (t^cc) i s plotted with a solid line. The parameters used to model the SN jet are 
E sn = 10 51 - 5 erg, tj = 10 s, t v = 10" 1 s, r 6 = 10 5 and e e = e B = 0.1. 



5.1. Meson suppression factors 

High-energy pions, kaons and muons produced by pnj and pp interactions do not all 
decay to neutrinos as electromagnetic (synchrotron radiation and IC scattering) and 
hadronic (irp and Kp interactions) cooling mechanisms reduce their energy. Muons 
are severely suppressed by electromagnetic energy losses and do not contribute 
much to high-energy neutrino production. Suppression factors for pion and kaon 
decay neutrinos are discussed next. 

The synchrotron and IC cooling times may be combined into a single electro- 
magnetic cooling rate as ig" 1 = + ijxj • F° r IC cooling in the Thomson regime 
t' lc w tf and in the KN regime t' lc 3> t' syn . The electromagnetic cooling time scales 
for mesons may be estimated assuming t' em w t' syn for simplicity as 

„ f 0.002 (E'JGeV)-i ( ^,0.5^,1^,-1 \ , , 

^-lO.S^/GeV)- 1 * { sb^E^ ) U4j 

The hadronic energy losses for mesons is similar to the proton energy losses by pp 
interactions in Eq. (|12|l with the same irp and Kp cross-section of 3 x 10 -26 cm 2 . 
The corresponding hadronic cooling time scales for mesons are 

c d - io- 5 - 4 ( lu ; J ^' '1 *, (15) 



E, 



51.5 



with AE' = Q.8E'. 

The electromagnetic (t'^. K cm ) and hadronic (t' had ) cooling time scales for mesons 
along with their decay times boosted by the respective Lorentz factors, t^. K doc , are 



February 5, 2008 7:55 WSPC/INSTRUCTION FILE snc 



10 Razzaque, Meszdros, & Waxman 




10-7L- — — — - — — ^ — 

10 2 10 3 10 4 10 5 

Meson energy E' K (GeV) 

Fig. 3. The Hadronic (i^ad)' electromagnetic (i^-ifem) ant ^ total (t' .jr coo {) cooling time scales 
for mesons: it and K as functions of comoving energy in a SN jet model by dashed, dotted and 
dot-dashed lines respectively. The meson decay time scales boosted by the respective Lorentz 
factors (t' .if dec ) are also plotted with solid lines. The parameters used to model the SN jet are 
E srL = 10 51 - 5 erg, tj = 10 s, t v = 10" 1 s, F b = 10 5 and e e = e B = 0.1. 



plotted in Fig. The total cooling time scale t'^. K cool = l/(C ; i, cm + C;i,had) is 
first dominated by the hadronic and then by the electromagnetic cooling channel. 
The ratio t'^. K coo \/t'^. K dec determines the suppression of mesons before they decay 
to neutrinos. 

From the condition t' cm — £ had , one may roughly define a break energy as 

= f 423/e B ,_i GeV (pion) 
br \ 6.5 x 10 4 /e B ,-i GeV (kaon), 1 ' 

above and below which the mesons cool by electromagnetic and hadronic interac- 
tions respectively. The corresponding suppression factor may be defined from the 
ratios t' em /t' dec and t' had /t' dec as 



W-^(E'jE' n>hl )-Pe B ,^ _ f 1 ; E'^ K < E' n . KM 
lO- 2 - 57 ^/^)-^,-! > p -\2 ] E% rtK > E' n . K>hl . 



(17) 



Another break energy may be defined from the condition t' dec = t' had below which 
mesons decay to neutrinos without any suppression. However, this energy: 24 and 
177 GeV respectively for pions and kaons in the comoving frame is low. The observed 
energy on Earth in both cases is below the detection threshold energy as will be 
discussed shortly. 



I 
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5.2. Neutrino fluence from individual sources 

If the shock accelerated protons could travel unimpeded from the SN jet at a lumi- 
nosity distance Dr,, then their isotropic equivalent fluence on Earth would be 

T p ob = g£|(l±f) (is) 

27rD i i? pi0b ln[£ , p :0 b,max/^p,ob,min] 

Here -E p . b = E p /(1 + z) and £j i0 b = tj(l + z) are the energy and time related in 
the observer's and local rest frames for the source location at redshift z. Of course, 
defelction by magnetic fields as well as other interactions along the way prevent 
such direct arrivale of the protons. 

Nearly all protons are expected to convert into n + (ir^ ) by jry (pp) interactions 
in the SN jet because r pj ;§> 1 (r pp 3> 1) in Eq. lfTU|l and assuming the charged 
pion multiplicity is ~ 1 from pp interactions as a conservative estimate. With a 
100% branching ratio for pion decay to neutrinos, one may define a multiplicative 
factor f n ~ 1 for protons which will produce neutrinos via pions. A similar factor 
fx ~ 0.063 may be defined for kaons as a product of multiplicity of ~ 0.1 
from pp interactions and 63% branching ratio for kaon decay to neutrinos. Protons 
produce neutral pions (kaons) and charged pions (kaons) with equal probability 
in both the %r) and pp interactions. For simplicity one may also assume that a 

carries 1/4 of the charged pion (kaon) energy: E v ~ E w> k/4, from roughly 
equipartition of energy between the final decay products of 7r ± and K . From these 
considerations one may estimate the neutrino (of one flavor) fluence on Earth as 

8 2irD L E u ob ln [E u ,ob,iaax/ E Ut ob,mm] \ E u ,Tr;K,bi ) 

per SN burst using Eqs. (|T7|) and (JTHJ). Here (3 = 2 (1) for E„. oh (l + z) greater 
(less) than E v ^;K,bT = ^bE' v . KhT /A from Eq. (JTBJl. The pre- factors in Eq. ifTYJ) 
are represented by the parameters r/^ = 10 -125 and r\K — 10~ 257 respectively for 
pions and kaons. For a typical ice Cherenkov detector such as IceCube, the threshold 
neutrino detection energy is E^.th ~ 10 E u %s GeV. The neutrino energy range is 
then ^, ob ,min - ^,ob,ma X = £„,th - 0.05r b £ pmax /(l + z) » 10 25 - 10 5 5 /(1 + z) 
GeV in the observer's frame on Earth. 

As mentioned earlier, 7r + 's are produced by pp or pj interactions by shock ac- 
celerated protons of all energies. In this case ir + decay fluence from a SN at a 
distance ~ 20 Mpc (D^ « 10 25,8 cm, z ~ 0), e.g. from the Virgo cluster, would be 

^,ob~10 ^ 10 2. 5 ) Dl M -i ' 

10 25 <£:^ob/GeV< 10 5 - 5 , (20) 

"Denoting the quantities in the observer's frame with subscript "ob". 
f Kaons are produced by pp interactions only. 
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from Eq. (|19fl with f v = 1. The neutrino break energy from pion decay: -Ev.7r.br ~ 
E u ,th ~ 10 2 ' 5 GeV in this case. A similar expression may be derived for neutrino flu- 
ence from kaon decays. However, pp interactions are overwhelmed by pj interactions 
in the energy range E' p ps 10 4 -10 5 2 GeV (Fig.EJ. Hence kaon and the correspond- 
ing decay neutrino production is expected to be suppressed in the energy range 
E v ps 10 3 9 -10 5 1 GeV. The neutrino break energy from Eq. (Q2J is E v>K ,br ps 10 47 
GeV and the fluence is 

■f* „ m -U.7( h ^ oh \ ^51.51 6,0.5 R v _i 2 . 

R _ J 1 ; io 2 - 5 < ^.ob/Gcv < io 3 - 9 

P - \ 2 ; 10 51 < E v , oh /GeV < IO 5 5 , ( > 
on Earth from kaon decays in a SN jet at a distance ~ 20 Mpc. 



5.3. Diffuse neutrino flux 



The diffuse neutrino flux is calculated by summing over fluences from all slow-jet 
endowed SNe distributed over cosmological distances in Hubble time. The SNe rate 
follows closely the star formation rate (SFR) which can be modeled, as a function 
of redshift per unit comoving volume!^ as 



0.32h 70 exp(3.4z) 



M Q yr _1 Mpc" 



(22) 



exp(3.8z) + 45 

70/i7o km s" 1 Mpc -1 is the Hubble constant. For a Friedmann- 



Here Hq 

Robertson- Walker universe, the comoving volume element is 



dV_ 

dz 



AirD 2 c 



1 + z 

and the relation between z and the cosmic time t is 



(23) 



(dt/dz)- 1 = -H (l + z)y/(l + n m z)(l + z) 2 - n A (2z + z 2 ). (24) 

For the standard ACDM cosmology, tt m = 0.3 and = 0.7. 

The number of SNe per unit star forming mass (/ sn ) depends on the initial mass 
function and the threshold for stellar mass to produce SN (M* ~ 8 M Q ). A Salpeter 
model 4>(M) oc M~ a with different power-law indices can generate different values 
for / sn , e.g. « 0.0122 M Q ' for a = 2.35 and ps 8 x IO" 3 Mq 1 for a = 1.35.E2EDI 
We adopt the model in Ref. EH which corresponds to the local type II SNe rate 
n sn (z = 0) = f sn p*(z = 0) ps 1.2 x 10~ 4 /i 3 yr~ 4 Mpc~ 3 agreeing with dataPH 

The distribution of SNe per unit cosmic time t and solid angle covered on the 
sky can be written as 

d 2 N sn _ n sn (z) dV _ h sn (z)D 2 L c 
dtdtt ~ 4tt dz ~ (1 + z) 2 



(25) 



I 



I 
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Fig. 4. Diffuse muon neutrino (vp) flux models from core collapse SNe assuming all of them 
Km, = 1 in Eq. EJ] are endowed with slow jets as discussed here. The pion decay neutrinos are 
produced via pp and p-y interactions over the whole range of shock accelerated proton energy. The 
kaon decay neutrinos are produced via pp interactions only and may be suppressed (thick dashed 
line) where p7 interactions are dominant over pp interactions. The thin dashed line corresponds to 
the kaon decay neutrino flux if pp is dominant over the whole proton energy range. The parameters 
used to model the SN jet are the same as in Figs.|5]and|!|] Diffuse fluxes arc unlikely to be detected 
because of the huge atmospheric background. 



We assume a fraction £ sn < 1 of all SNe involve jets and a fraction l/2r| of all such 
jets are pointing towards us. The observed diffuse SNe neutrino flux, using Eqs. 
(|P3j) and is then 



diff 



2Tl 



dz 



d 2 N sn 
dtdtt 



$~ ^,ob (-^i/,ob ) 



£snfTT,Kr]Tr,K cEj 



32tt 



P 2 

-'Vob ^0 



dz 



n sn (z)/(l 



hi [E v 

X (E^th < -Ev.ob < -Ev,ob,: 



E v>oh [l + z] 



E, 



(26) 



with (3 specified as before. 

The diffuse flux from all cosmological slow-jet SNe is plotted in Fig. 21 by 
numerically integrating Eq. (|26[) . assuming the maximal fraction £ sn = 1. The pion 
(kaon) decay flux is plotted with thick solid (dashed) line(s). The exponential sup- 
pression in the kaon decay flux curves is due to the lack of kaon production 
by pp interactions, as pry interactions, which do not produce kaons, are dominant 
in the particular proton energy range corresponding to this neutrino energy range, 
-Ev,ob(l + z) k 10 39 -10 51 GeV. The light dashed curve corresponds to the kaon 
decay flux if pp interaction dominates over the whole proton energy range for 
comparison. The kaon decay neutrino fluxes are first smaller and then larger than 
the pion decay neutrino flux. This is because of the kaon's heavier mass and shorter 
decay time compared to pions (see Fig. |3J). 
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The atmospheric v^, flux from pion and kaon decays (conventional fl ux) 
compatible with AMANDA datd^l is plotted with the following parametrizatiorP"! 



<T) i O.012£;- O 2 b 74 /(1 + 0.002^ ob ) ; E v , oh < 10 5 8 GeV 

' S 3 - 8 ^ob 7 /(l + 0.002£„, ob ) ; E v , oh > 10 5 - 8 GeV. 



Also shown are the cosmic ray bounds (WB limits) on the diffuse neutrino flux^l 
It is unlikely that neutrino detectors can measure the SNe diffuse fluxes plotted in 
Figure 01 as they are below the atmospheric background. However, individual SNe 
in nearby galaxies may be detectable, as discussed in the next section. 

6. Neutrino Detection Prospects 

There are approximately 4000 galaxies known within 20 Mpc distance. At the stan- 
dard rate of I SNu = 10~ 2 /yr/10 10 blue solar luminosity for average galaxies, the 
estimated SN rate is > 1/yr. The SN rate in the starburst galaxies, such as M82 
and NGC253 (at distances 3.2 and 2.5 Mpc in the Northern and Southern sky, 
respec tive ly) is ~ 0.1/yr, much larger than in the Milkyway or in the Magellanic 
cloudsPSl Very strong neutrino signals in future kilometer scale neutrino detectors 
are expected from these nearby SNe, over a negligible atmospheric background, 
using temporal and positional coincidences with optical detections. 

The directional sensitivity of the Cherenkov detectors is best for muon neutrinos, 
which create muons by charge current neutrino-nucleon (yN — > fiX) interactions. 
Muons carry ~ 80% of the incident neutrino energy. It emits Cherenkov light as it 
travels faster than the speed of light in the detection media. Photo multiplier tubes 
(PMTs) buried in the medium (ice, e.g. in the case of IceCube) can detect muons 
by gathering their Cherenkov light. The effective detection area of a Cherenkov 
detector depends on the arrival direction of the neutrino and the energy of the 
muon. In principle it can be larger than the geometrical area (A^ct ~ km 2 in the 
case of IceCube e.g.) as muons can be produced outside the instrumented volume 
and travel inside. To achieve a good pointing resolution a muon should hit at least 4 
PMTs, for reconstructing its track unamiguously, strung on different vertical strings. 
The PMT spacing is 17 m vertically and 125 m horizontally. With a muon energy 
loss of ~ 0.2 GeV/m and a PMT efficiency of ~ 30%, the threshold energy for muon 
detection is (3 x 125 m) x (0.2 GeV/m)/0.3 « 250 GeV as a conservative estimate. 
This corresponds to a neutrino threshold energy E v ^ w 10 2 ' 5 GeV as used before. 
The IceCube detector is expected to have an angular sensitivity of 1° for muon 
tracks of energy < 1 TeV coming from a zenith angle < 140° and it gets better at 
higher energyP3l 

The number of muon events from a nearby jetted SN may be calculated as 

N„ = A dct f " m ' lX V{E V , e)F*„. K (E v )dE u , (28) 

where the neutrino fluences are given in Eqs. H20|) and (|2ip. The full detection 
probability V(E V ,&) depends on the source's angular position (9) and hence on the 
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Fig. 5. Cumulative muon events from flux models from a jetted core collapse SN in M82 and 
NGC253 (lower and upper sets of lines) as functions of observed muon energy in IceCube. The pion 
and kaon decay neutrino flux models are described in Eqs. H2UI . 1211 and in the caption of Fig, El 
(also in the main text). The neutrino events from pion decay flux models are larger (smaller) than 
the kaon decay flux models at lower (higher) energy. This is because more pions are produced than 
kaons but pion decay to neutrino is suppressed due to pion energy loss at higher energy compared 
to kaons with heavier mass and shorter decay time (Fig.0. 



Earth's shadowing effect, as well as the energy dependent vN cross section!""' The 
cumulative number of muon events from M82 and NGC253 are plotted in Fig. [S] 
per individual SN with its jet pointing towards Earth . The flux models are from 
pion (solid lines) and kaon (dashed lines) decays in the SN jet as in Fig.|H The two 
kaon decay models are due to the protons in the SN jet producing kaons in their 
whole (thick dashed lines) and partial (thin dashed lines) energy range (see Fig.[2J). 
The lower and upper sets of lines at — 10 3 GeV correspond to the SN in M82 
and NGC253 respectively. The timing uncertainty in the optical detection of SN 
is *~ 1 day. The corresponding atmospheric background events within 1° angular 
resolution is 0.07 for both M82 and NGC253 in the same energy range as in Fig. El 

Neutrino fluxes of all three flavors {y^, u e , v T ) on Earth should be equal because 
of oscillations over astrophysical distances. However, only v^'s are emitted from the 
sources under consideration and the total number of neutrino events (including v e 
and v T ) will remain the same as the events calculated here. The lack of good 
directional sensitivity for the v e and v T events may prevent obtaining their positional 
coincidence with the SN. Timing coincidences of v e and v T events with events, 
however, may still be useful to verify the neutrino oscillations at these energies and 
test their common origin. 

At the quoted rate for M82 and NGC253, a SN from one of these galaxies 
would be expected within five years. However, only 1/5 of them would be visible 
in neutrinos due to the 1/2T% beaming effect of the jet. IceCube may detect one 
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such SN in 25 years of its operation. The real situation is not as pessimistic as for 
the brightest neutrino events expected from M82 and NGC253. Other nearby spiral 
galaxies (M31, M74, M51, M101, etc, and the Virgo cluster) will also have SNe. 
For a hypothetical SN at 20 Mpc with its jet pointing towards Earth, the number 
of neutrino events is « 15/km 2 from pion and kaon decay fluxes combined. With 
the suggested rate of > 1 SN/yr within 20 Mpc, IceCube may detect > 3 muon 
events/yr after beaming correction from the jetted SNe. 

7. Conclusions 

While a core collapse SN in a typical galaxy is a rare event, their rate of occurence 
within 20 Mpc could be more than 1/yr, and the physics and astrophysics one can 
learn from such an explosion is enormous. A buried slow jet from the collapsing 
core of the supernova progenitor star is an attractive possibility for solving the long 
standing problem of how to understand the ejection of the envelope in SN explosions, 
by re-energizing the shock wave through energy deposition by the jet. Alternatively 
it could be that only a fraction of core collapses leads to such jets. These hypothetical 
mildly relativistic jets in SNe may be related to the ultra relativistic jets thought 
to be responsible for long duration GRBs, which are thought to originate from 
the core collapse of massive progenitor stars, some of which have been positively 
associated with observed envelope ejection supernova events. While all typical core 
collapses should produce 10 MeV thermal neutrinos, the presence of a jet would 
allow proton acceleration by shocks, and produce 1 TeV non-thermal neutrinos. 
Detection of these high energy neutrinos by upcoming Cherenkov detectors would 
be a smoking gun signal of a SN jet, and would allow one to study the conditions 
inside a collapsing star which may not be possible otherwise. 

Acknowledgments 

Work supported by an NSF grant AST0307376 and NASA grant NAG5-13286. 
References 

1. For a recent review see: M. Delia Valle, |arXiv:astro-ph/05u45T71 

2. A. MacFadyen and S. E. Woosley, Astrophys. J. 524, 262 (1999) 
[arXiv:astro-ph/9810274| . 

3. E. Ber ger, S. R. Kulkarni, D. A. Frail and A. M. Soderberg, Astrophys. J. 599, 408 
(2003) arXiv:astro-ph/0307228 . 

4. J. S. Bloom, S. R. Kulkarni, F. A. Harrison, T. Prince, E . S. Phinney and D. A. Frail, 
Astrophys. J. 506, L105 (1998) [arXiv:astro^p h/9807050"l . 

5. B. M. S. Hansen, Astrophys. J. 512, L117 (1999). 

6. A. I. Macfadyen, S. E. Woosley and A. Heger, Astrophys. J. 550, 410 (2001). 

7. E. Waxman and A . Loeb, Phys. Rev. Lett. 87, 071101 (2001) 
|arXiv:astro-ph/0102317| . 

8. E. Ra mirez-Ruiz, A. Celotti and M. J. Rees, Mon. Not. Roy. Astron. Soc. 337, 1349 
(2002) arXiv:astro-ph/0205108 . 



February 5, 2008 7:55 WSPC/INSTRUCTION FILE snc 



High energy neutrinos from a slow jet model of core collapse supernovae 17 

9. J. Granot and A. Loeb, Astrophys. J. 593, L81 (2003) |arXiv:astro^p h/0305379]. 

10. K. Nomoto, P. A. Mazzali, T. Nakamura, K. Iwamoto, I. J. Danziger, F. Patat, in 
Supernovae and gamma-ray bursts: the greatest explosions since the Big Bang, May, 
1999, eds. M. Livio, N. Panagia and K. Sahu. STScI symposium series, Vol. 13, pg 
144 (2001). 

11. J. Granot and E. Ramirez-Ruiz, Astrophys. J. 609, L9 (2004) 
arXiv:astro-ph/0403421 . 

12. L. Wang, D. A. Howell, P. Hoflich, and J. C. Wheeler, Astrophys. J. 550, 1030 (2001). 

13. D. C. Leonard, A. V. Filippenko , A. J. Barth and T. Matheson Astrophys. J. 536, 
239 (2000). 

14. For a recent review see: A. Mezzacappa, arXiv:astro-ph/0010580 

15. A. M. Khokhlov, P. A. Hoflich, E. S. Oran, J. C. Wheeler, L. Wang, and 
A. Yu. Chtchelkanova, Astrophys. J. 524, L107 (1999). 

16. B. Zhang and P. M eszaros, Int. J. Mod. Phys. A 19, 2385 (2004) 
|arXiv:astro-ph/0311321| . 

17. E. Waxman and J. N. Bahcall, Phys. Rev. Lett. 78, 2292 (1997) 
|arXiv:astrc-ph/9701231| . 

18. P. Meszaros and E. Waxman, Phys. Rev. Lett. 87, 171102 (2001) 
|arXiv:astro-ph/0103275| . 

19. S. Razzaque, P. Meszar os and E. Waxman, Phys. Rev. D 68, 083001 (2003) 
|arXiv:astro-ph/0303505| . 

20. S. Razzaque, P. Meszaros and E. Waxman, Phys. Rev. Lett. 93, 181101 (2004) 
[Erratum-ibid. 94, 109903 (2005) ] |arXiv:astro- ph/0407064"l . 

21. S. Ando and J. F. Beacom, arXiv:astro-ph/0502521 

22. K. Hirata et al. [KAMIOKANDE-II Collabora tion] , Phys. Rev. Lett. 5 8, 1490 (1987). 

23. F. Cei, Int. J. Mod. Phys. A 17, 1765 (2002) arXiv:hep-ex/0202043 . 

24. J. Ahrens et al. [IceC ube Collaboration], Astropart. Phys. 20, 507 (2004) 
|arXiv:astro-ph/0305196] . 

25. I. Sokalski [ANTARES Collaboration], arXiv:hep-ex/0501003 

26. C. L. Fryer and R. Dupuis, in 30 Years of Discovery: Gamma- Ray Burst Symposium. 
Edited by E. E. Fenimore and M. Galassi. AIP Conference Proceedings, Vol. 727, 371 
(2004). 

27. W. Zhang, S. E. Woosl ey and A. I. MacFadyen, Astrophys. J. 586, 356 (2003) 
|arXiv:astro-ph/0207436| . 

28. S. Razzaque, P. Meszaros and E. Waxman, Phys. Rev. Lett. 90, 241103 (200.3) 
|arXiv:astro-ph/0212536] . 

29. C. Porciani and P. Madau, Astrophys. J. 548, 522 (2001) arXiv:astro-ph/0008294 . 

30. L. He rnquist and V. S pringel, Mon. Not. Roy. Astron. Soc. 341, 1253 (2003) 
|arXiv:astro-ph/0209183| . 

31. P. Madau, M. Delia Valle and N. Panagia, Mon. Not. R. Astron. Soc. 297, 17L (1998) 
|arXiv:astrc^p h/9803284 . 

32. J. Ahrens et al. [AMA NDA Collaboration], Phys. Rev. D 66, 012005 (2002) 
|arXiv:astro-ph/0205109| . 

33. M. Thunman, P. Gon dolo and G. Ingelman, Astropart. Phys. 5, 309 (1996) 
arXiv:hep-ph/9505417 . 

34. A. Alonso-Herrero et al, Aston. J. 125, 1210A (2003). 

35. see e.g. S. Razzaque, P. M eszaros and E. Waxman, Phys. Rev. D 69, 023001 (2004) 
|arXiv:astro-ph/0308239| . 



I 



I 



